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Every d x al bipartite system is shown to have a large family of undistillable states with nonpositive 
partial transpose (NPPT) . This family subsumes the family of conjectured NPPT bound entangled 
Werner states. In particular, all one-copy undistillable NPPT Werner states are shown to be bound 
entangled. 
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Introduction: Maximally entangled bipartite states 
form an indispensable resource in several quantum in- 
formation processing situations 0, 0, EJ ■ In reality, 
however, readily available states may be mixed and less 
than maximally entangled due to a variety of reasons in- 
cluding influence of the environment. Thus distillation, 
the process by which some copies of (nearly) maximally 
entangled pure states are extracted from several copies 
of partially entangled mixed states using local quantum 
operations and classical communication (LOCC), is of 
singular importance. 

Rapid progress has been achieved in this regard, be- 
ginning with the works of Popescu [j| , Bennett et al. , 
Deutsch et al. Q , and Gisin U|jJ . An important question 
that presented itself at an early stage of this development 
was this: Can every entangled state be distilled? The 
Horodecki family answered this question affirmatively in 
a significant particular case: Every inseparable state of 
a pair of qubits can, give n sufficiently many copies, be 
distilled into a singlet [lfj . 

Subsequently they proved a result that applies to arbi- 
trary to x n bipartite systems 0] : A state that does not 
violate the Peres-Horodecki PPT (positive partial trans- 
pose) criterion [l2L Hij | can not be distilled. Since in- 
separable PPT states can not be distilled, they are said 
to possess bound (i.e., undistillable) entanglement. The 
first examples of such undistillably entangled PPT states 
were constructed by Horodecki , and several families 
of PPT bound-entangled states have been presented since 
then. The role PPT bound-entanglement could play as a 
quantum resource has also been studied. 
The Problem of NPPT Bound Entanglement: 
Given a bipartite state p, let (p TB )® n = (p® n ) TB denote 
the partial transpose of the state of 77 identical copies. 
Then the necessary and sufficient condition for distilla- 
bility of p is that the inequality 0] 
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be violated by a Schmidt rank two vector in the n- 
copy Hilbert space, for some 77. Thus in order to be dis- 
tillable, the state should be NPPT (i.e., the state should 
have non-positive partial transpose) . What has remained 
an open problem is this: Are there undistillable NPPT 
states? Does undistillability imply PPt[T^? 



Clearly, violation of the above inequality for a partic- 
ular 77 = no implies violation for all 77 > 770. Given an 
NPPT state p, if it satisfies this inequality for all Schmidt 
rank two states, for a particular 77, we say that the state 
is 77-copy undistillable or, equivalently, that the partial 
transpose p Ts is 77-copy 2-positive. Thus p is undistill- 
able if p Ts is 77-copy 2-positive for all 77. 

While Horodecki et al. pointed out in a subsequent 
work 0| that NPPT bound entangled states, if they exist 
in nature, should be found among the one-copy undistil- 
lable NPPT Werner states ^3] > it m &y be fair to say that 
this problem attained its present state of fame only when 
two leading groups produced independently, and about 
the same time, analytical and numerical evidence |l8l Hif 
for its existence in the context of the Werner family of 
states [More recent numerical attempt has been under- 
taken in Ref. [HI]- As the analytic part of the evidence 
it was shown, for every 77, that there is a corresponding 
range of parameter values over which the NPPT Werner 
state remains provably 77-copy undistillable. But the 
range itself becomes rapidly smaller with increasing 77, 
and was not proved to remain nonzero as 77 — * 00. It 
should be added that no one-copy undistillable NPPT 
Werner state was shown to be 77-copy distillable either. 
Similar evidence in the multipartite case has been con- 
sidered by explicit construction of 77-copy undistillable 
NPPT states, for every 77 > 1 [2lj|. 

It should be noted for completeness that some systems 
are known not to support NPPT bound entangled states 
at all. These include 2 x 77 dimensional systems 0, 0] 
and all bipartite Gaussian states jl^l • 

Shor et al. have used the conjectured existence 
of NPPT bound entanglement to prove nonadditivity 
of bipartite distillable entanglement. And Eggeling et 
al. |2J| have related the existence of NPPT bound entan- 
glement to the connection between the sets of separable 
superoperators and PPT-preserving channels. Further, 
Vollbrecht and Wolf |2fij have shown that an additional 
resource in the form of infinitesimal amount of PPT 
bound-entanglement can render any one-copy undistill- 
able NPPT state one-copy distillable. 

In a more recent paper Watrous has constructed, 
for every 77 > 1, states which are 77-copy undistillable, and 
yet are distillable. This is a surprising and important 
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result, for it shows that n-copy undistillability, even if 
n is very large, does not by itself prove undistillability. 
The burden this finding places on numerical evidence for 
NPPT bound entanglement is evident. 

We conclude this brief summary of the present status 
of the conjectured existence of NPPT bound entangle- 
ment by noting that the conjecture itself seems to enjoy 
the confidence of researchers in quantum information the- 
ory, even though the evidence presented so far has been 
assessed differently by different authors mum 

In this Letter we prove that any d x d bipartite sys- 
tem with d > 3 has a fairly large family of NPPT states 
which are undistillable. As will be seen, this family is 
much larger than the Werner family of conjectured NPPT 
bound entangled states, but it turns out to be convenient 
to begin our proof with the Werner family. 
Proof of Existence of NPPT Bound Entangle- 
ment: We may define the one-parameter family of 
Werner state p a in d x d dimensions through the par- 
tial transpose of p a : 

p T B = Id _ daP ( 2 ) 

Here P is the projection on the standard maximally en- 
tangled state: 

1 d 

P=|*><*|, |*> = -^= (3) 

These states as defined are not normalized to unit trace, 
but this does in no way affect our considerations below. 

Nonnegativity of p a forces on a the restriction — 1 < 
a < 1. This allowed range for a divides into three inter- 
esting regions [ljl [l9j : 

- 1 < a < PPT, separable, 
d 

- < a < 1, NPPT, one-copy distillablc, 

~ < oi < -, NPPT, one-copy undistillable. (4) 

Clearly, it is the last mentioned range which is of interest 
for the issue on hand. In this range p a is NPPT, yet p a TB 
is 2-positive. Since p a being NPPT implies that p® n is 
NPPT, for all n, the issue really is whether (p a TB )® n is 
2-positive as well. We show that it indeed is, for all n. 

Let us denote by S^ 1 ' the collection of all pure and 
mixed states of Schmidt rank (number) < 2 2^ . This is a 
convex set whose extremals are all pure states of Schmidt 
rank 1 or 2. Let 9 — (6\, 62, • •■ , O4) be a d- tuple of angles, 
and consider the subgroup of U(d), the ci-dimensional 
unitary group, consisting of diagonal matrices Ug\ 

(U e ) M = 5 kl e ie *. (5) 

This is the standard maximal abelian subgroup Ua of 
U(d). Now average each element a of over the local 



group U A ®U* A : 

g e -> a = J dOUe® Ug a U\ <x> Uj . (6) 

Here dO = d9\d9i ■ ■ ■d9d 1 with 9 k running over the in- 
terval [0, 2tt], independently for every k. This group- 
averagingprocess is analogous to the diagonal twirl op- 
eration [18|. The resulting images of elements of con- 
stitute a new convex set f^ 1 ), where the superscript on 
S and fi reflects our intention to extend these consider- 
ations to several copies of the Werner state. 

The reduction in complexity achieved by going from 
to r^ 1 ' should be appreciated. For instance, whereas 
has a CP^ 1 x CP^ 1 worth of product states 
among its extremals, the product states among the ex- 
tremals of f^ 1 ) are precisely d 2 in number; these are the 
standard computational basis states \k, I). Similarly the 
pure states of Schmidt rank = 2 among the extremals of 
fid) necessarily have one of the forms ai|ll) -I- 012 1 22), 
#2 1 22) + /3 3 |33), and 7i|ll) + 7 3 |33); each of these three 
sets is a Bloch-sphere in worth. And, therefore the maxi- 
mally entangled rank-two states are precisely three circles 
worth in number, being of the form -^(|11) + e lSl \22)), 

^(|22) +e"»|33)) > or ^(|33) + e«»|ll)). 

Now 2-positivity of p^ B is equivalent to the demand 
that Tt(p Tb a) > 0, for all a S In view of the 

Ua ® U\ symmetry of p Tn , this is equivalent to the de- 
mand that Tr(p TB a) > 0, for all a G fi (1) and m view 
of the convexity of f^ 1 ), it is sufficient for the extremals 
of Q^ 1 ' to meet this demand. This convexity argment is 
in essence a recognition of the fact that the real-valued 
expression Tr(p^ B sigma) is linear in a. 

We thus find that the minimum of Tr(p Ts cr) over fK 1 ) 
equals (1 — 2a), showing that p Tn is 2-positive for all 
a < |. The minimum is achieved by -^=(|11) + |22)), 

^(|22) + |33», and ^(|33) + |11», and by no other 
states. 

We may note in passing that the local symmetry of 
p Ts is not U A <8> U%, but the full group U{d) ® U{d)*. We 
wish to go beyond the Werner states later in this Letter, 
and for this reason we have based our analysis on the 
subgroup Ua <8> U A , rather than on the full group. 

We now move on to consider p a ® p a , two copies of 
the Werner state p a . The set is constructed as 
the collection of (d 2 x d 2 )-dimensional states of Schmidt 
number < 2. From the convex set we obtain 

fi( 2 ) by averaging each a G over the local group 

(Ua ® U%) ® (Ua ®U* A ). It is to be understood that 
the first (Ua ® U A ) factor acts on the Hilbert space of 
the first copy and the second on that of the second copy, 
independently. 

The extremals of fl^ are readily enumerated. The 
rank-1 extremals of fi' 2 ' are necessarily of the form 
|fc, l)i ® \i,j)2- They are (i 4 in number. The rank two 
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states are of two types. 

Type I: (a\kk)i + <8> \ij) 2 , 

|fc/)i® («|m>2 + /3|jj) 2 ); 

Type II : a\kk) x <8> + 0\ll)i ® |jj>2- (7) 

At the risk of sounding repetitive, we emphasize that this 
is a complete enumeration of the extremals of Q( 2 \ 

As with the single copy case 2-positivity of p a (g> p a is 
equivalent to the demand Tr[(p^ B <g> p„ b )<t] > 0, for all 
a € f^ 2 ), which in turn is equivalent to the demand that 
this condition be met by all the extremal states of $}( 2 \ 
The type-I rank-2 states being products across the two 
copies, cannot bring out genuinely two-copy properties if 
any, and thus we are left with only the type-II states to 
examine: 

{^ii\pI b ®pl B \*n) 

= \a\ 2 (kk\p^ B \kk)(ii\p^ B \ii) 

+m 2 (u\Pi B \u)(jj\pi B \jj) 

+a*(3{kk\p T ?\ll){ii\p T a B \jj) 
+af3*{ll\pZ°\kk){jj\pl°\ii). (8) 

Now, 2-positivity of p Tn is equivalent to the Schwartz 
inequality 

{kk\pl-\ll)<[{kk\pl B \kk){ll\pl-\ll)]K (9) 
for rank-1 states. Use of this inequality in Eq.(8) proves 
<i&Il\ P l B ®pl B \*Il)>0. (10) 

That is one-copy 2-positivity of the (U A <8> U A ) invariant 
operator p Ts implies its two-copy 2-positivity. And we 
have proved 

Theorem 1: All one-copy undistillable Werner states, 
that is all p a 's in the entire range g < a < |, are two- 
copy undistillable. 

To move on to the n-copy case, assume that p^ B is 
(n - l)-copy 2-positivc. That is (ip(pl B )® (n " 1) |^) > 
for all rank-2 states of the (n — l)-copy Hilbcrt space. We 

wish to prove that this implies (p^ 13 )^^ is 2-positive. 

To this end, form fi(") by averaging each a £ over 
the local group (Ua <S> UX)® n , with one (Ua ® U A ) factor 
acting on the Hilbert space of each copy independently. 
Again, the extremals of f^ n ) consist of rank-1 and rank-2 
pure states. The rank-1 states are of the form <8> 
K 2 , :? 2 ) <8> ■■■■\i n ,jn}- That is, these are tensor products of 
computational basis states, one picked from each copy, 
and thus are d 2n in number. 

The rank-2 extremal (or pure) states in fi( n ) are of two 
types, as in the two-copy case: 

Type I: = ® \i, j), 

Type II : \V H ) = |^) ® + \<t> 2 ) ® \j,j), (11) 



where \ip) is a (n — l)-copy rank-2 state, and \4>i), 1^2) 
are (n — l)-copy rank-1 states. Since type-I states have 
a product structure across the copies, nonnegative ex- 
pectation values of (p^ B )®" in respect of type-I rank-2 
states follows directly from the assumed (n — l)-copy 2- 
positivity, the additional copy offering nothing new in 
the type-I case. The same (n — l)-copy 2-positivity is 
equivalent to the validity of the Schwartz inequality for 
(n — l)-copy rank-1 states, and this in turn implies the 
nonnegativity of the expectation values of (Pa B )® n f° r 
type-II states. We have thus proved, by induction, 
Theorem 2: The one-copy 2-positive (pJ B )®"'s in the 
entire parameter range \ — a — \ are n-copy 2-positive, 
for all n. That is, these one-copy undistillable NPPT 
Werner states p a are n-copy undistillable, for all n, and 
hence are bound entangled. 

Finally, it should be evident that the only property 
of p a , apart from its NPPT and one-copy 2-positivity 
properties, used in our analysis is its (Ua®Ua) invariancc 
or, equivalcntly, the (Ua ® U A ) invariance of its partial 
transpose. It follows that our conclusions apply to all 
states with these properties. That is, 
Theorem 3: Every one-copy undistillable NPPT state 
in d x d dimensions is bound entangled if it possesses 
(Ua®U a ) symmetry. 

This is the main result of this Letter. It shows that 
the family of NPPT bound entangled states in d x d di- 
mensions, for any d > 3, is much larger than what might 
have been anticipated. This point is worth illustrating. 

It is easily seen that in 3 x 3 dimensions the most gen- 
eral (Ua g" U A ) invariant P Tb has the form 
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Clearly, p has to be positive semidefmite in order to be 
a valid density matrix. This demand is equivalent to the 
conditions: (i) all the diagonal elements pij > 0, and 

PijPji > \ z ij\ f° r a U i 3- The NPPT requirement 
demands that the 3x3 submatrix 

Pll — 212 —23! 
~Z*2 P22 -2 23 
_ -231 -^23 P33 

should be nonpositive, and the 2-positivity demand is 
equivalent to the three inequalities P11P22 > 1 2:12 1 2 , 
P22P33 > I223I 2 , and P33P11 > |2 3 i| 2 - The NPPT demand 
and the 2-positivity demands thus involve only the six 
parameters pn, P22, P33 and z\ 2 , Z23, 231. The phases of 
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the complex z-parameters can be tuned by (local) change 
of phases of the basis vectors on the A and B sides (to 
be precise it is sufficient to carry out the changes on one 
side only), but the argument of z 12 Z23 £31 is invariant 
under such gauge transformations. 

Thus our family of NPPT bound entangled states in 
3x3 involves, when normalized to unit trace, 12 param- 
eters; eight coming from the diagonals pij, three coming 
from the magnitudes of the z-parameters, and one gauge- 
invariant phase. These are canonical parameters, and do 
not take into consideration parameters arising from local 
unitary transformations. 
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